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Abstract 

We analyse a minimal supersymmetric standard model (MSSM) taking a minimal flavour vio- 
lation (MFV) structure at the GUT scale. We evaluate the parameters at the electroweak scale 
taking into account the full flavour structure in the evolution of the renormalization group equa- 
tions. We concentrate mainly on the decay B s — > and its correlations with other observables 
like BR(fc — ► 57), BR(6 — > sl + l~), AMb s and the anomalous magnetic moment of the muon. We 
restrict our analysis to the regions in parameter space consistent with the dark matter constraints. 
We find that the BR(i? s — > can exceed the current experimental limit in the regions of 
parameter space which are allowed by all other constraints thus providing an additional bound on 
supersymmetric parameters. This holds even in the constrained MSSM. Assuming an hypothetical 
measurement of BK(B S —> ~ 1CP 7 we analyse the predicted MSSM spectrum and flavour 
violating decay modes of supersymmetric particles which are found to be small. 
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I. INTRODUCTION 



The Minimal Supersymmetric Standard Model (MSSM) is the simplest supersymmetric 
(SUSY) structure that includes the SM gauge group and matter content. However, this def- 
inition does not unambiguously fix a single supersymmetric model. As is well-known, SUSY 
cannot be an exact symmetry of nature and must be broken. Therefore, to specify completely 
the model, it is necessary to fix the soft breaking terms: this amounts to 124 parameters 
at the electroweak scale. Although these parameters would be fixed in the presence of a 
truly fundamental theory, from the point of view of the effective theory they are only con- 
strained by low energy experimental data. In fact, our freedom in this enormous parameter 
space is already severely limited by phenomenological constraints. The accumulating con- 
cordance between the SM expectations and the vast range of experimental results in FCNC 
and CP violation point toward a low energy new physics which is, at least approximately, 
flavour blind jjj. In this spirit, the so-called Minimal Flavour Violating (MFV) MSSM is an 
MSSM where the only non-trivial flavour structures present are the usual Yukawa couplings 
while all sfermion masses and trilinear couplings are completely family universal Q . Notice 
that realistic flavour models have been constructed where the deviations from this minimal 
flavour structure are rather small 3|. 

Several interesting flavour changing processes can take place for large values of tan/? 
at a rate much larger than the corresponding expectations from the SM. A particularly 
interesting process in this scenario is the decay B s — > which has been shown to 

receive enhancements of even three orders of magnitude with respect to the standard model 
expectation, for ,a rg e tan , and s m a„ psendoscalar H lggS .asses |fc . This p r oee 8 s has b ee„ 
carefully analysed in previous studies. In fact, two different points of view have been taken 
in these analyses. On the one hand, several groups examined this decay in the framework of 
an effective MFV MSSM model defined at the electroweak scale where masses and mixing 
an g ,es are unrelated flflflflHHQ. In these worfcs ,ar g e ennancen-ent factors are 
relatively easy to find because several correlations with other observables are lost. On the 
other hand, various other works repeated this analysis in a MFV MSSM defined at the GUT 
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scale in terms of a reduced number of parameters 
Most of these analysis take place in the framework of a constrained MSSM (CMSSM) defined 
at Mqut in terms of mo, mi/2, Aq, tan f3 and sgn(/x) and neglecting flavour mixing in the RGE 
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evolution. In principle it is possible to consider different family-universal masses for different 
representations under the gauge group. However, to our knowledge, in the literature so far 



only t 

mass 



re masses of the Higgs multiplets are allowed to be different from a common sfermion 



23, 



241 . |25| with the only exception of Reference |26[ . In any of these GUT scenarios, 
several relations between the sfermion masses, gaugino masses, Higgs masses and mixing 
angles are obtained and although large enhancements of the B s — > branching ratio 

are still allowed, they are clearly more difficult to obtain. 

In this work we analyse this decay in the framework of a completely generic MFV MSSM 
defined at the GUT scale. In this general MFV scenario, the soft masses and trilinear 
terms for different representations under the SM gauge group are different, although family- 
universal. This scenario can be realised in different ways, as RGE-induced splittings between 
different multiplets in the running from Mpianck to Mgut in 5 77( 5) or flipped SU (5) models 
(27 1 . or in models of direct unification at the Planck-scale Another difference with 

respect to previous analysis is the fact that, for the first time, we use the full two loop renor- 
malization group equations (RGE) [3] with the complete flavour structure of the different 
flavour matrices to determine the Yukawa couplings and soft SUSY breaking parameters at 
the electro-weak scale. Furthermore, we take into account all the relevant constraints in our 



analysis including the dark matter relic abundance |30| together with the updated results 



on processes like B s -B s mixing 
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321 ] . muon anomalous magnetic moment |33| ^6 — > sl + l~ 
|34ll35j and the LEP and Tevatron constraints from searches of SUSY particles All the 
analysis is done with updated values of the input parameters and in particular we use the 
last value of the top quark mass . Fi nally , we take into account the latest constraints on 
the BR(7? S — > from D0 and CDF 38j] . In this framework, we study the neutral Higgs 

effects in the decay B s — > and the B s mass difference, comparing the results of the 

CMSSM with our generic MFV scenario. We take special care in checking the compatibility 
of these processes with other tan /3 enhanced observables like the b — > 37 decay or the muon 
anomalous magnetic moment. 

The structure of the paper is as follows. We start with a discussion of the origin of Higgs- 
mediated flavour changing neutral currents (FCNC) and study their phenomenological effects 
in the B s — > decay and the B s mass difference which are enhanced by additional tan /3 

factors. In section IIHI we analyse the correlations of these Higgs-mediated processes with 



other tan/3 enhanced observables like the b — > 57 decay, the muon anomalous magnetic 
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FIG. 1: Flavour changing coupling of the up Higgs-boson, H u , to the down-type quarks. 

moment and the b — > sl + l~ decay. In section IIVI we describe the procedure used in our 
numerical analysis and present our results for the B s —>■ fi + n~ decay in the various scenarios. 
Furthermore we discuss the various correlations with other observables. Section 0is devoted 
to the discussion of the expected collider phenomenology if the decay B s — > is found 

with a branching ratio close to the present bounds. Finally in section IVII we present our 
conclusions. 



II. HIGGS MEDIATED FCNC'S 

The MSSM is classically a type-II two-Higgs doublet model. This means that one of 
the Higgs fields, H u , couples only to the up quarks while the other, H d , couples only to 
the down quarks. In this way, dangerous tree-level Higgs-mediated FCNC's are avoided. 
However, the presence of the /i-term in the superpotential breaks the symmetry protecting 
this structure. Therefore quantum corrections modify this type-II structure and generate 
couplings of H\ to up quarks and of Hi to down quarks, hence reintroducing Higgs-mediated 
FCNC's. We use in our numerical analysis the full expressions for the one-loop vertex and 
self-energy corrections as given in Reference Q (see also Q)- All our numerical results 
are valid for any value of tan (5. However, in the following we summarize the main features 
of Higgs-mediated FCNCs, relevant for our analysis, in the limit of large tan (3. A more 
detailed discussion can be found in Reference . In Fig. d we present an example of 
FCNC coupling of the up- type Higgs, H u , to the down quarks through a stop-chargino 
loop. There are also similar couplings mediated by gluino loops and the corresponding 
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couplings to up quarks. The effective Yukawa Lagrangian at 1 loop is, 

C Yuk = -d^Yid^Hi + d u (AY d ) .. d K] Hl + u hi Y%u Rj H u + u u (AY% u R .H* d , (1) 

where Y$ and Y$ include also the 1 loop corrections to the i^-down-quark and if u -up-quark 
couplings ^(|. After spontaneous breaking of the electroweak symmetry, this Lagrangian 
gives rise to new contributions to the quark mass matrices. Working in the basis where the 
tree-level down-quark Yukawa couplings are diagonal and neglecting the subdominant loop 
corrections in the large tan/3 limit, the contributions to the down quark mass matrix are: 

= tJ| du (Yfa + tan/3 (AY d ) J d Kj + H.c, (2) 

where we can see the appearance of a tan/3 enhanced correction to the down quark mass 
due to the presence of the wrong-type Higgs vacuum expectation value (vev). In the same 
basis, we obtain for the neutral Higgs couplings to the down quarks [lo| : 

C s° = ~du ( n m 4 5 iJ + ( AYd ) ■ ■ 4) d ^ S ° + Hc -> ( 3 ) 



x v d (I + 9 tan/?) d 13 V Jij ~ u ) 

with S° = (H ,h ,A°,G ), xf = (cos a, — sin a, i sin/?, —i cos (3) and xf = 
(sin a, cos a, —i cos /3, —i sin (3). Here, Wij denotes the running quark mass that can be ex- 
tracted by experiments, the parameter ej accounts for the tan f3 enhanced corrections to the 
down quark eigenvalues 10, 3^ and it is given by 



€j = e + \K 3j \ 2 y 2 e y , 

/ 2 2 s 

e = -2a s — ^— H : ' 



— 1 o ' 2 1' 

d7T rrig \ m~ m~ J 

l07T /i \ fi Z fi Z ) 

where the loop function is defined in Appendix D. Using the fact that 1 -£^2(1, 1) = — 1/2 and 
taking the limit m- g ~ // ~ m^ a ~ — A t , we obtain eo — sign(/x) -0.012 and £3 ~ sign(/i) -0.015. 
Finally ( AY ') is given by: 

\ Jij 

(AY d ) . . = jfc, (eotfy + y t 2 # 3 * K 3j e Y ) , (5) 



1 We have that H2(x,x) goes from —0.18 for x — 4 to —1.7 for x = 0.1 



5 



with y d . being the eigenvalues of the tree-level Yukawa matrix. 

After diagonalizing the full one-loop down-quark mass matrix, the couplings to the neutral 
Higgs-bosons are given by 



C 



s° 



-d 



(4 + tjxf) Sij + 



U Vv d (1 + e J -tan/3) ^ d ' " v ' v d (1 + l, tan/5) 2 

£y ( 4 - zf tan/5) W 5° + H.c 



(6) 



with A 



o 



and Aq = K^K^j. Hence, from this equation we can see that we do have Higgs- 



mediated FCNC's after the diagonalization of the one-loop mass matrix and in particular 
tan j3 enhanced off-diagonal couplings proportional to xf which come from the rotation di- 
agonalizing the one loop Yukawa couplings. In the following we study the phenomenological 
consequences of these FCNC Higgs couplings in the down-quark sector 2 . 



A. The decay B s — > fi + fi 

The present experimental bound on this decay from the CDF and D0 
tions is, 

BR(B S -> fj+fM-) < 1.0 x KT 7 95%C.L. 



38, 



collabora- 



(7) 



This process is specially sensitive to the presence of Higgs mediated FCNC's. Using the 
Lagrangian given in Eq. @ we obtain the dominant tan /^-enhanced neutral-Higgs terms to 
this decay which contribute mainly through two effective operators: 

2G F 



n 



a 



cff 



-K1K, 



C s m b (6 R s L ) (/ I J + C P m b (&rSl) (l 75 I) 



V2 2irsm 2 6 w " tb ts 

The dominant contributions to the Wilson coefficients Cs and Cp in the large tan/3 limit 



are 



Cj 



z A tan 3 /3ff 2 (m 2 L //i 2 ,m 2 R //T 



4M^ fi (l + e tan/3)(l + e 3 tan/3) 
m M m 2 A tan 3 (3 H 2 (m\J\x 2 , m\J 
AMw fi (l + e tan/3)(l + e 3 tan/3) 



sin(a — (3) cos a cos(a — j3) sin a 



M 2 H0 



Ml 



(9) 



2 Other papers study the large tan effects in r decays |4JJ |41| |42| , flavour changing decays of the Higgs 
bosons 0> and collider processes 0, 0, . 
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For large tan (3 we can neglect the M 2 contribution and, taking M B o ~ Ma, we find that" 
Cs — Cp. Therefore, the leading contribution to the 6rSl — > amplitude is given by: 

m 6 m M K* b K ts A t tan 3 (3 



A(b R s L -> fi ff 



oc m. 



H 2 (m\J ii 2 , m\J y) (10) 



* (1 + e tan/?)(l + e 3 tan/5) /i 
Using this formula, we obtain an estimate of the branching ratio in the large tan f3 region 

as, 



BR [B s -> ii+ff 



2.2 x 10" 



1.5 ps 



230 MeV 



m t 



175 GeV 



tan p 
50 



350 GeV 



A| (H 2 {m 2 Jfi 2 ,m 2 Jff 
n 2 (1 + e tan/5) 2 (l + e 3 tan/3) ; 



(11) 



and we see that the branching ratio for this decay scales like tan 6 /3/m\. As we will show 
below, given the new experimental bounds from CDF and D0, this process is already able 
to probe a sizable part of the large tan f3 parameter space. 



B. B® s mass differences 



The CDF and D0 collaborations [32| have recently presented the first experimental results 
on B,, mass-difference: 



AM Bs = (l7.33l^(stat) ± 0.07(syst)) ps" 1 , 
while the mass difference is 

AM Bd = (o.500+°;^( s tat) ± 0.07(syst)) ps" 1 . 



The basic formula for the B q mass difference is 



□ 



AM 



B a 



G 2 F M 2 W 
6n 



M BqVB F B B Bq \K tq \ 2 \F t V 



with 



n 



S (m 2 /M 2 w ) + ^C n v e L w L (/i) 



1 



+ ^CV VRR (//) + A LR ^1 LR (^) + ^P^V) + 



pf LL [cf L (/i) + cf RR (ii)} + p 2 sll 



(12) 



(13) 



(14) 



Cf LL (//) + Cf RR (/i)l . (15) 



However, a difference between these masses is possibie for light pseudoscalar Higgs and can give rise to 
sizable effects 0, Ql • In our numerical calculation, we use the full expression and we have found that 
for Ma > 250 GeV, Cs — Cp is a very good approximation. 
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The different Wilson coefficients are defined in and r = 0.985 describes the 0(a s ) 

radiative corrections to & (m,/tfw) in the SM. 

As shown in Refs. [10, 50], the same flavour-changing Higgs couplings generating the 
leading contribution to the decay B s — > fi + fi~ at large tan/? induce a Higgs-mediated 
double penguin contribution to B sd -B s ^ mixing. In the large tan/? regime, the dominant 
contributions to the B S; d mass differences are given by the usual SM box and by the double 
penguin contribution 



AM Bq oc \K tq \ 



q(m 2 /M 2 s _ m b m q mt A t 2 tan 4 /? 

J a\"H Ivl w) ' 2 nj2 7T~i I a\2n , ~ I ml 2 2 x 
7r Mw (1 + e tanp) (1 + e 3 tan/?) fJ- Tn-A 

H 2 (m 2 t ' J ii 2 , m\ R //i° N N " 



(16) 



with q = s,d and S (x) is defined in the appendices. From this expression we can see that 
the double penguin contribution is suppressed by an additional factor of Gf with respect to 
the SM contribution due to the fact that it is a two loop contribution. However, it has a 
strong dependence on tan /? which can lead to a huge enhancement compensating this loop 
suppression. Moreover, these double penguin contributions have always opposite sign with 
respect to the Standard Model box contribution and, thus reduce the mass difference for 
large tan /?. This effect is specially important in the case of AMg s yielding an additional 
constraint in the large tan/? regime. Notice that in the case of AM# d , the double penguin 
contribution is the same as in AMb s with an additional suppression of md/m s . Given that 
this contribution competes with the same SM function So(x t ) it is clear that this process 
will not be relevant to constrain the SUSY contributions. 



III. CORRELATIONS WITH OTHER PROCESSES 

Given the strong dependence of B s — > n + n~ on rriA and tan f3, we can expect that 
other processes, specially tan j3 dependent processes, receive large SUSY contributions if 
BR(5 S — > is at the 10~ 7 level. A typical example is the decay B —>■ X s ^ whose 

measured BR agrees very well with the SM prediction. In turn this provides a stringent 
constraint to SUSY contributions and given that the chargino amplitude is proportional 
to tan/? this holds in particular for the large tan/? regime. The decay B — > X s £ + £~ is 
closely related and the combination of both processes has been recently shown to eliminate 
the possibility of changing the sign of the dipole Wilson coefficient C-j in MFV scenarios 
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5jJ,|52|. Another important constraint is provided by the muon anomalous magnetic moment 
which is also tan/? dependent and constitutes a very important constraint specially for 
fi < 0. 



A. B - X sl 



The Standard Model contributions to the decay B — ► X s 7 are known at NLO |53j,|54j and 
the excellent agreement with the experimental results constrain strongly any extension of 
he SM. Charged-Higgs contributions in two-Higgs doublet models are also known at NLO 



53 



55 



561 ] . In the MSSM a complete NLO calculation is still missing but the most important 



contributions have been calculated in Refs. 
process mostly through the dipole operators 



57 



58, 



59 



6fl. The MSSM contributes to this 



-"Dipole 



V2 



K tb K ts 



C 7 (ji) ■ ^s L a, v b R F^ + CM ■ ^s La T^a, u b R pG a ^ 



16tt 2 



The experimental world average from the CLEO 
collaborations is given by 



Belle 



62, 



(17) 

and BaBar 0, 



BR(B -> X s7 )^ >1 . 6GeV = (3.55 ± 0.24°j° 9 10 ± 0.03) x 10~ 4 
This has to be compared with the standard model prediction 0] 



BR[B Xtf] = (3.61 



+0.24 
-0.40 



mc/mb 



± 0.02 CKM ± 0.24 



par am 



± 0.14 scale ) x 10- 



(18) 



(19) 



In the numerical analysis we utilize the formula presented in Ref. [67] in which the branching 
ratio is explicit given in terms of arbitrary complex C-j and C§. 

The two Wilson coefficients, Cj(fi) and Cg(/i), receive contributions from W-boson, 
charged Higgs and chargino diagrams. It is very interesting to discuss the relative signs 
of these three contributions. Charged Higgs and W-boson diagrams have the same sign and 
interfere always constructively. In contrast, the chargino contribution can have either sign 
depending on the sign of the /i parameter. In the large tan /3 region the relative sign of the 
chargino mediated diagram is given by — sign(A t/ u). In the scenarios studied in this paper 
sign(A t ) = — sign(Mi/ 2 ) as can be seen in table IV of 2^|. For this reason, /i > implies 
destructive interference while /i < implies constructive interference of the chargino and 
W-boson contributions. 
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Given that the SM contribution essentially saturates the experimental bound, \x < 
scenarios are tightly constrained and the sum of charged Higgs and chargino contributions 
is restricted to be a small correction. However, for positive \x the chargino contribution has 
opposite sign with respect to the W-boson and charged Higgs ones. In this case, a cancelation 
between sizable charged Higgs and chargino contributions allows to keep the total Wilson 
coefficient within a small margin from the SM expectation. In principle it is possible for 
the chargino to cancel both the charged Higgs and W-boson contributions and generate a 
total Wilson coefficient of equal modulus and opposite sign to the usual SM contribution. 
This is consistent with experimental data because the BR(5 — > X s 7) depends essentially 
on |C| (/^b)| 2 and therefore is not sensitive to the sign of these coefficients. However, as we 
will see later, combining this process with the decay B — ► X s £ + £~ it is possible to constrain 
the sign of the CV Wilson coefficient 51 



and exclude this possibility in the generic MFV 
scenarios studied. 

The complete expressions for the LO contributions to the Cj and C% contributions that 
are used in our numerical evaluation are well-known and can be found in Ref. 0,0]- For 
the following discussion, it is enough to examine the Cj coefficient as Cs behaves in a similar 
way, 

C?W) = £f~/i (m 2 t /m 2 w ) , (20) 

Ct(m w ) = I {1 ^ ^ A (mllm^) + f 2 fa/nfe) • 

( 1 (eo + e 3 )tan/? , 2 t Y t' Y tan 2 /? \\ 

\ l + e 3 tan/? ^ (1 + e 3 tan/?) (1 + e tan/?) J J ' 1 ' 

6 a,a'i=Ia=l n tb n ts m j( ± < Z V 



m Xl 1 



Q{^',a)I* H {a,a)I ^ (m^Jml^J | , (22) 



m w a/2cos/?(1 + e 3 tan/3) 
where, for simplicity we have only included in this equation the leading contributions in 



the limit of unbroken SU(2) x U(l) symmetry 



10]. 



In the numerical analysis we use the 



complete expressions with full flavour dependence. The parameters e' and e' Y are given by 



The main tan j3 effects are then present in the charged Higgs couplings and in the chargino 
couplings to right-handed down squarks. K aq G^ a ' kS>I represents the coupling of the chargino 
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/ and the squark k to the left-handed down quark q; and m q /(-\/2mw cos (3)K aq H( a > k > 1 the 
coupling of the chargino i and of the squark k to the right-handed down quark a. T hese 
couplings, in terms of the standard mixing matrices defined in the appendices are (70| 



\ y2m w sin p J 

H (a,k)i = -U I2 T k v a L . (24) 

The explicit expressions for the loop functions can be found in the appendices. 

In the large tan/? regime, it is interesting to expand the Cj coefficient assuming the 
smallness of the off-diagonal entries with respect to the corresponding diagonal elements both 
in the stop and chargino (more precisely in M^M~) mass matrices. In this approximation 
we obtain, 

7 Qm^ sin 2 3 1 + e 3 tan (3 m 2 ± m 2 ^ — m~ R 

^ fiA t m 2 tan/? 5 m 2 ^ 

~ Qm 2 w sin 2 3 1 + e 3 tan f3 12 m\ ^ ' 

where in the second line we replaced the combination of loop functions by its value in the 

limit in which stop left, stop right and chargino have the same "average" SUSY mass, m|. 

In this limit the SUSY contribution to the amplitude is 

4Gj? r „,ff r „ ff * 5 uA t m 2 t tan 3 

A(b R - Sl7 ) (x -£k*K« m b ^ ' * / . (26) 
V2 72 m| (1 + e 3 tanp) 

If we compare now this amplitude with the corresponding amplitude of the B s — > fi + fi~ 
decay, Eq. ([26)1 . we see that, although both them are enhanced by powers of tan/3, these 
two amplitudes behave differently: ^4(6rSl — > scales as tan 3 j3/m 2 A while the b — > s^y 

amplitude scales as tan/?/m 2 (assuming /x, A t and mj of the same order). Taking into 
account this difference we can analyze if it is possible to obtain a large contribution for 
B s — > fi + fi~ while at the same time satisfy the stringent bounds on the b — > 57 decay. 
Obviously, the answer to this question depends on the sign of \i. 

In the B — > decay with \i < chargino contributions to C7 and Cg have the same 
sign as charged Higgs and W-boson contributions. If we require a sizable neutral Higgs 
contribution to B s — > n + ji~ we need a light pseudoscalar Higgs, which implies also a light 
charged Higgs, and large tan/3. A light charged Higgs enhances while large tan/3 

enhances C7 . As both contributions compete to fill the narrow space left over by the SM 
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contribution in the B — > X s 7 amplitude these two requirements are clearly incompatible. 
The only possible solution is to increase the squark masses (in particular the stop mass) 
to suppress the chargino contribution to CV while keeping as low and tan/3 as large 
as possible. In this case we are essentially in a two-Higgs doublet model without SUSY 
and from the existing constraints on the charged Higgs mass 3, 56] we infer that the 
pseudoscalar must be equally heavy. Moreover, the constraints from the muon anomalous 
magnetic moment, that are analyzed below, make the /i < situation even more difficult. 
Therefore it is not possible to reach very large values of BR(i? s — > as we will see in 

the numerical analysis of the next section. 

The case fi > is much more interesting phenomeno logically. A large BK(B S — > /i + /i~) 
still requires a light charged Higgs and large tan/3, however the chargino contributions 
interfere now destructively with the charged Higgs and W-boson ones and thus these two 
conditions are fully compatible because |C? + Cf + | < 0.1CV SM . As we will see in the next 
section, this cancelation must be quite precise and this constraint remains very strong. In 
the flavour-blind models defined at the GUT scale, in particular in the CMSSM, is not 
always easy to obtain the required values of m^, and tan/? to achieve this cancelation. 
The extent to which this cancelation can be realized will be addressed in the next section 
where we perform a full analysis of the parameter space. 

B. B — > X s l + l~ 

The effective Hamiltonian that describes b — > sit transitions in the Standard Model is 

AC ( 10 6 \ 

H eS = 4-Kt b K* ts £ Qin) O t (ji) + £ C iQ {n) O lQ ^) + C b OM , (27) 

V2 \i=l | _i=3 / 

where we adopt the same definitions as in Ref. |52j ]. The operators OiQ and Ob are required 
only for the inclusion of electro-weak corrections. In particular the most relevant Wilson 
coefficients are C 2 ([ib) (that does not receive sizable new physics contributions), CV(/i;,), 
C 8 (/ife), C 9 (/ib) and Ci (/ife). The explicit definitions of 7 and 8 are given in Eq. (fTTj) and 
the semileptonic operators read 

9 = ^P (s Llfl b L ) {hn) , O w = ^ (Sl-tM (hV) ■ (28) 

47T 47T 

The branching ratio BR(B -> X S U) is known at NNLO in QCD and at NLO in QED (only 
terms enhanced by large logarithms - i.e. \ogm w /m b and logrrib/mi - are included). 
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Intermediate cc resonances produce peaks in the dilepton invariant mass (4m 2 < s < m 2 B ) 
spectrum and disrupt quark-hadron duality. For this reason one has to consider the low- 
and high-s regions separately. The very low-s region (s < 1 GeV 2 ) is dominated by the 
quasi-real photon emission and contains the same amount of information we already extract 
from b — > 57. In the high-s region (s > 14.3 GeV 2 ) the integrated branching ratio is quite 
small; moreover, the l/mj, expansion breaks down and reliable predictions for the spectrum 
suffer large uncertainties. The low-s region (1 GeV 2 < s < 6 GeV 2 ) is dominated by the 
Wilson coefficients C 9 and Ciq, has sizable integrated branching ratio and is very sensitive 
to the C 7 — Cg interference. In the following, we utilize the latter to put constraints on the 
Wilson coefficients £7,8,9,10- 

We calculate the integrated branching ratio in the low-s region following Ref. [5: 



BR/ 



2.2306 - 0.0005 1(R 10 ) + 0.0005 I(R 10 R* 8 ) + 0.0906 X{R 7 ) + 0.0223 X{R 7 R* 8 ) 
+0.0050 1(R 7 R* Q ) + 0.0086 I(R 8 ) + 0.0258 I(R 8 R* 9 ) - 0.0615 1(R 9 ) 
-0.5389 K(R W ) + 0.1140 K(R 7 ) + 0.0154 n(R 7 R* w ) + 0.0687 K(R 7 R* 8 ) 
-0.8414 K(R 7 R* 9 ) - 0.0036 K{R 8 ) + 0.0019 n(R 8 R* w ) - 0.0980 TZ(R 8 R* 9 ) 
+2.6260 U(R 9 ) - 0.1074 TZ(R 9 R* 10 ) + 10.6042 l^o] 2 + 0.2837 |i2 7 | 2 



+0.0039 |i? 8 | 2 + 1.4645 \R 9 



x 10~ 7 , (29) 



where Ri = C\ ot (/i ) / 'Cf M (fi ) are the next-to-leading order Wilson coefficients (i.e. they do 
not include 0(a s ) corrections) at the high scale normalized to their SM values. When dealing 
with the dipole operators we utilize the standard definition of the scheme-independent effec- 
tive coefficients [52]. In the numerical analysis we set to zero supersymmetric contributions 
to the NLO matching conditions. 

The standard model and supersymmetric contributions to the leading order matching 



conditions of the semileptonic operators are 52j (we do not give explicit formulae for the 



gluino, neutralino and neutral Higgses contributions because their effects are negligible com- 
pared to the charged Higgs and chargino ones): 



= ^Y(x t ) + W(x t ) + \ - t log ^ , (30) 

CFo = -i-Y{xt) , (31) 
s w 

H ± 4Sjy — 1 X t . 2 . 2 . 1 

85^ t&n 2 (3 J5K 11 H ! 18 tan 2 /? 



Or = I2 I L / 5 K7^ ± ) + -.TZZT-a h^Vm 2 H± ) , (32) 
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c 



10 



9 



^10 



1 



f 5 (m 2 /m 2 H ±) 



Ssfy tan 2 /3 

6 2 j^eff 7>-eff * 

ly ab ly a's Q(a',A)I*Q(B,a)J 



(33) 



<*,<*' A,£=l Z,J=1 ^Hfe zv ts 

c 2 (m 2 ± , m? A , m 2 A ) (r^ L r^ L ) ^ 5/j + 



4s 2 



<? 2 



^f^c (ml A , m 2 f , mj^CMn 



- c 2 (m 2 A , m 2 ± , m 2 ± ) <5ab Vj\ V n 



l^fr(ml f /ml A )S AB 6 u 



+—m 2 v d 2 (m 2 ± , m 2 ± , m? , m? )5 AB V,*!^ 



1 



(34) 



2 T^eff T^cff * 
y y y n al) Ji a's ^(o ! ',A)/*^(( J B,a)J 



o 2 ^cff^cfT* 

c 2 (m| ± , m 2 A , m 2 A ) (IV L 4 L ) ab 5 U + 2 



+ 



2m 2 w d 2 (m 2 ± , m 2 ± , m? A , m 2 x ) - c 2 (m? , , m 2 ± , m 2 ± ) 



(35) 



A scan of the CMSSM and MFV parameter space results in tiny deviations of the Wilson 
coefficients C 9 and C 10 from their SM values: |C 9 SUSY /C 9 SM | < 0.3/1.65 and |Cf ^sy^sm | < 
0.8/4.45. 

This decay has been observed by Belle js^J and BaBar j^j]. In the low-s region the 
experimental results read as 



BR(5 
BR(73 



XjTt 
XJ + t 



(1.493 ± 0.504+°;^}) x 10" 6 (Belle) 
(1.8 ± 0.7 ± 0.5) x 10~ 6 (BaBar) . 



(36) 
(37) 



This leads to a world average 



BR(E -> X s 



1.60 ±0.51) x 10" 



(38) 



That is in perfect agreement with the SM prediction, that reads: 



BR(5 -> X & 



0+1 



Ism 



1.58 ± 0.07 scale ± 0.06 mt ± 0.025 mc ± 0.015 



nib 



±0.023 Qs(A / z ) ± 0.015 CKM ± 0.025 BRsi 
(1.58 ±0.10) x 10" 6 . 



x 10" 



(39) 



In order to assess the impact of this measurement on the supersymmetric parameter space, 
we perform a model independent analysis. As a first step we extract the allowed ranges for 
the Wilson coefficients C-j and Cg from the measured B —>■ A s 7 branching fraction. This 
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C 7 < C 7 > 




-6 -4 -2 2 -8 -6 -4 -2 



Rg Rg 

FIG. 2: Impact of the BR(i? — > X s l£) measurement on the Rg-Rio plane. The black disk is the 
region accessible in an MFV model, i.e. it includes the possible SUSY contributions both in the 
CMSSM as in the generic MFV scenario that we analyse here. 

gives two thin stripes in the C-j — C$ plane that correspond to the Cj < and CV > 
scenarios. We further restrict Cg to be within a factor of 10 around the SM value. As a 
second step, we consider the two strips separately, vary C-j and Cg in the allowed range, and 
plot the resulting constraint in the Cg — Cio plane. The whole analysis is performed at the 
90% C.L. and theory errors are included. Our results are summarized in fig. The most 
striking feature of these plots is that, in the C-j > scenario, the SM point is excluded and 
contributions in MFV supersymmetric scenarios are not large enough to bring theory and 
experiment into agreement; thus we conclude that in such models the sign of CV has to be 
SM-like ( i.e. negative). For a similar discussion see also jod) ]. 



C. Muon anomalous magnetic moment 



In SUSY theories. a„+ receives contributions via vertex diagrams with X°~P> and x —v 



loops 111 I . The chargino diagram strongly dominates in almost all the parameter space. For 
simplicity, we will present here only the dominant part of the chargino contribution (the 
complete expressions that we use in the numerical simulation can be found in Ref. |3]): 



mi 
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where the loop function ^3 is given in the appendices. In fact, if we expand this function 
in the small off-diagonal elements of the chargino matrix we obtain for the dominant tan /3 
enhanced contribution, 

V- 327r 2 m 2 m 2 > l 4i J 

where we have chosen the convention M2 > as we do in the rest of the paper. The most 
relevant feature of Eqs. (pfUj) and (j4*Tj) is that the sign of 5a*+ is fixed by sign[Re(L r i 2 Vii)] = 
— sign [Re (/i)]. The experimental result for this observable is at present 72], 

a* xp = 11 659 208 (6) x 1(T 10 , (42) 

while the theoretical expectations within the SM are 



a 



SM. 



e+e") = 11 659 181 (8) x 1(T 10 , (43) 
G s M ( r ) = 11659196(7) x 1(T 10 . 

Comparison with the experimental results implies that these results strongly favour the 
fi > region in a MFV scenario. This is specially true for the theoretical prediction based 
on e + e~ annihilation that is smaller than the experimental resul t by 2.7a. In the case of 



resul t by 
1.4a [TJ 



the prediction based on r decay the difference is reduced to 1.4a [72j but it still requires a 
positive correction and disfavours strongly a sizable negative contribution. Therefore this 
has important consequences in all our observables. As we saw in the previous section, \i > 
is precisely the phenomenologically interesting region for the combination of the B s — > fi + fi~ 
and B — > Xgj decays. On the other hand the a M measurement makes even more difficult to 
find interesting effects in B s — > [i + for /x < 0. In the numerical analysis, presented in the 
next section, we use a conservative a M bound at 3a with the e + e~ based estimate of the SM 
contribution and we still find that with \x < we can only find BR(5 S — > /i + /i~) at the level 
of 10~ 8 with very heavy SUSY spectrum. In the /i > region, the conservative 3a bound 
allows for a vanishing supersymmetric contribution. 



IV. NUMERICAL RESULTS 



As said in the introduction, in this work we concentrate on a version of the MSSM with 
MFV soft-breaking terms at the GUT scale. This means that the only non-trivial flavour 
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structures in the model are the usual SM Yukawa couplings; all soft-breaking masses are 
family universal and the trilinear couplings are proportional to the corresponding Yukawa 
matrices. The RGE evolution between the GUT and electroweak scales introduces small 
non-universal entries in the soft-breaking terms. 

In this framework we analyse two different versions of the MSSM: the CMSSM and 
the most general MFV MSSM. The full spectrum of the CMSSM is determined (assuming 
vanishing SUSY phases) by five parameters: My 2 , m , A , tan/3 and sign(/x). Therefore, 
the SUSY parameters at the electroweak scale, including couplings, masses and mixing 
angles, show several interesting correlations. Furthermore, these five parameters are strongly 
constrained by the experimental bounds on masses and different FCNC processes. In this 
restricted framework it is very interesting to investigate how large BR(i? s — ► oc 
tan 6 (3/m\ can be, while remaining in agreement with all the other experimental constraints. 

We can ask the same question in a more general flavour-blind MSSM, where we allow for 
different, but generation- diagonal, parameters for the fields with different gauge quantum 
numbers. Here the number of parameters is thirteen: Mi/ 2 , ttiq^ m u R fl, m d R ,oi m Z,0' m e fl ,o, 
m H d ,o, m H u ,o, Ai,o, Ad t o, Ae t Q, tan/3 and sign(/i), where the in the subscripts indicates that 
these parameters are defined at Mqut- The precise definition of Mqut and the numerical 
evaluation of the parameters at the electroweak scale will be described below. Clearly in 
this model we have more freedom than in the CMSSM, implying milder correlations between 
the masses of different particles; thus, we can expect larger rates for the Higgs-mediated 
processes consistent with other FCNC constraints. 

In the numerical analysis presented below we define our model, CMSSM or most general 
MFV, in terms of the GUT scale parameters listed above. We scan the values of these pa- 
rameters in the following ranges: M1/2 < 1 TeV, m^o < 2 TeV (i = Q, Ur, cIr, L, 6r, H u , Hj), 
A 2 U>0 < 3(m| +m| K)0 +m| u)0 ), A\ Q < 3{m^ Q +m 2 iRfi +m 2 Hdfi ), A 2 efl < 3(m 2 ifi +m 2 gRfi +m 2 Hdfi ); 
|// 1 is calculated from the requirement of correct electroweak symmetry breaking. The bound 
on the Aifi parameters is set to avoid charge and/or colour breaking minima. 

At the low scale, we impose the following constraints on each point: 



Lower bound on the light and pseudo-scalar Higgs masses 
The LEP constraints on the lightest chargino and sfermion masses [36] 
The LEP and Tevatron constrains on squarks and gluino masses 
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Agreement with the experimental data on low ener 



y p rocesses like BRTB — > X .</y) 



lergy p roc 



34, 



35 



Mm, 



61 



62, 



63, 



|6J,|65|; 

• Correct dark matter relic abundance from stable neutralinos 

We consider MSSM with conserved R-parity and thus, the lightest neutralino (that coincides 
with the lightest supersymmetric particle - LSP) provides an excellent dark matter candidate. 
In fact, the recent precise determination of the cold dark matter density of the universe is a 
very strong constraint on the SUSY parameter space. In most of the parameter space of the 
MSSM the neutralino abundance generates a too large density of cold dark matter. Only 
certain regions of the parameter space, in which neutralinos are efficiently annihilated, are 
in agreement with the experimental results. 

The most precise results on the cold dark matter density come from the WMAP experi- 
ment 

Q CDM h 2 = 0.1047ig;gg at 68% C.L. (44) 

This constraint selects very narrow strips in the MSSM parameter space jz^j]. On the other 
hand it has been shown 76j that these allowed regions are extremely sensitive to small 
differences in the computation of SUSY masses. This is specially true in the large tan/3 
region that we analyze in this paper, where differences in the stability of radiative symmetry 
breaking or in the Higgs masses are due to the different approximations used in the RGE 
programs. In view of all this, we prefer to take a conservative attitude with respect to the 
dark matter constraint and, allowing the possibility of other cold dark matter candidates, 
impose only the loose 99% C.L. upper bound, 

ficDM^ 2 < 0.13 (45) 

In this section we analyze the different parameters (masses, mixing matrices, . . . ) that 
enter the B s — > fi + fi~ amplitude, and the impact of the various constraints that we have 
discussed above. We start with the analysis of this process in the CMSSM and then generalize 
it to the most general MFV model at the GUT scale. 



4 This constraint assumes a standard thermal history of the universe and can be easily evaded in non 
standard cosmological models. 



18 





5.110 10~ 4 


pole 
m t 


172.9 




0.1057 


m b (m b ) 


4.25 


m T 


1.777 


m z 


91.1876 


m u (Q) 


3 • 10~ 3 


Gf 


1.1664 -lO" 5 


m d (Q) 


7 • 10~ 3 


1/a 


137.036 


m s (Q) 


0.12 


Aa Ld 


0.02769 


m c (m c ) 


1.2 


af s (m z ) 


0.1172 


S12 = 0.224, s 23 = 4.13 • 10" 2 , 


si 3 = 3.63 • 10~ 3 , 5 = 1.13 



TABLE I: Numerical values of the SM input. Masses are given in GeV, for the leptons and the 



t quark the pole masses, for the lighter quarks the MS masses either at the mass scale itself, for 
c, b, or, for u, d, s, at the scale Q = 2 GeV. Sij are the sines of the CKM mixing angle and 5 the 
CKM phase. 

A. Calculation of the parameters at the electroweak scale 



In the numerical calculations, we use the following procedure. The masses and mixing 
angles at the electroweak scale are calculated in an iterative way: First the gauge and 
Yukawa couplings are calculated from the input values given in Table at the electroweak 



scale taking into account the shift from the MS to the DR scheme as well as the effect of 
the SUSY thresholds. In case of the SUSY thresholds the complete flavour structure of the 
sfermions is taken into account in the calculation. Afterward these couplings are evolved 
to the GUT scale which is defined by the requirement g\ = g 2 where g\ = J5/3g' is the 



properly normalized U (1) coupling and g 2 is the SUl{2) coupling. We do not require at M GU t 
that #3 is equal to g± = g 2 but assume that this differences is accounted for by unknown 
thresholds due to very heavy particles with masses of the order Mqut- At Mgut the SUSY 
breaking boundary conditions are set. At this scale the trilinear couplings are obtained 
multiplying the A-parameters with the Yukawa matrices taking into account the complete 
flavour structure. Also in the RGE evolution, the complete flavour structure is taken into 

rj 

account and we evolve the parameters at the two- loop order |29j. This clearly induces off- 
diagonal elements in the sfermion mass parameters. The calculation of the SUSY masses is 



carried out in the DR scheme and the formulae for the one-loop masses are generalisations 
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M Bs 


TBs 




f R 1/2 
JBsD B s 


r(B s ) 


pLR 


jbLR 


pSLL 


pSLL 


5.3696 GeV 


1.461 ps 


0.55 


(294 ± 33) MeV 


1.461 ps 


-0.71 


0.9 


-0.37 


-0.72 



TABLE II: Numerical values for the B-physics parameters using the notation of ref. [ic| . 
lattice parameter has been taken from Ref. 



The 



of the ones given in Re: 



presented elsewhere 



s 



taking into account the complete flavour structure and will be 



In the case of the /i parameter and t 



he masses o: 



81 



82 



the neutral Higgs 



83|, |8J, |85| are added. Note 



bosons the two-loop corrections as given in Refs. 0, 
that in the one-loop contributions we use the full flavour structure, while for the two-loop 
part we use the approximation that only the sfermions of the third generation contribute. 
The obtained spectrum is then used to recalculate the thresholds to gauge and Yukawa 
couplings restarting the complete procedure. This is repeated until the relative difference of 
all couplings and masses between subsequent iterations is at most 1CT 5 . 

The calculation of the Yukawas, in particular Y^, is important for the calculation of low 
energy observables. In the calculation of the Yukawa couplings (see Eq. ((TJ) we take into 
account the thresholds of supersymmetric particles as given in 10]. Due to the RGE running, 
off-diagonal elements in the squark sector are induced which remain non-zero when going 
to the super-CKM basis. For this reason, not only the chargino contributions to the off- 
diagonal terms but also gluino and neutralino contributions are taken into account although 
the later ones are negligible. The gluino contribution amounts typically to a 10% correction 
of the chargino contribution, which is precisely the order of magnitude one would expect for 
the two-loop corrections. 

For the calculation of the low energy observables, the complete formulas are used, in- 
cluding all contributions stemming from H + , xt> Xj> 9-> h°, H°, A in case of B physics 
observables and the contributions due to y^, in case of a M . More precisely, in case of 
B s — > fi + fi~ we combine the formulas of jlfl| . jso] ] and However, we use only the lowest 
order contributions to C$, Cp, C' s , C' P , C w and C[ because the complete one-loop QCD 
corrections to these Wilson coefficients are not available in the literature. In the case of 
AMp q we have used the formulas of Refs. rJ 88]. In the case of b —>■ s'j we have used the 
brmulas given in j3] and we have added the leading SM QCD corrections as described in 



3 



. For the SUSY contributions to the anomalous magnetic moment of the muon we have 
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used the formulas given in [89(. The hadronic parameters we use in all these calculations 
are collected in Table |Hj In most of the cases we neglect the uncertainty in these param- 

~ 1 /2 

eters with the exception of /b s -B b ' that has a large impact in the B s mass difference. In 
the calculation of the B s mass difference we compare below the results obtained fixing this 
parameter to the value in Table ITTl and the results varying fB s B B ^ within the one-sigma 
range f Bs B B f2 = (0.295 ± 0.036) GeV Q- 



B. B s -> n+ii- in the CMSSM 

As we have seen in previous sections, tan/3 enhanced Higgs FCNC at low energies are 
inversely proportional to the pseudoscalar Higgs mass and are essentially generated via a 
stop-chargino loop. In the context of a CMSSM defined at Mqut it is natural to ask what 
are the possible values of these masses and couplings (m^, tan (3, . . . ) that can be reached 
consistently with all the phenomenological bounds. 

We will analyze this problem discussing separately the case of /i > and fi < 0. In fact 
sign(/i) is a key parameter in the analysis of the phenomenology of the large tan (3 limit for 
two main reasons. First, as well-known and explained in detail section ITTT1 the sign of the 
SUSY contributions to several low energy processes like b — > 57 and the muon anomalous 
magnetic moment are determined by the sign of the /1 parameter; therefore the stringent 
experimental constraints tend to favour the positive sign of \i specially in the large tan/3 
limit. In addition, as we have seen in sectioning the sign of the chargino one- loop corrections 
to the down-quark Yukawa matrices and in particular of the eo and ij parameters are also 
fixed by sign(/i). This implies that, in the basis of diagonal tree-level down Yukawa matrices, 
these one-loop corrections tend to reduce or increase these diagonal entries depending on 
the sign of /i and hence they give rise to a different phenomenology. In the \x > case, 
we have that both eo and ij are positive; therefore, Eqs. (JJJ)-(jnj) imply that the diagonal 
Higgs couplings to the down quarks are smaller than the naive expectation, fn^/vi. As a 
consequence, for instance, FCNC couplings are smaller and effects of down-quark Yukawa 
matrices in the RGE's are reduced. 

The requirement of correct electroweak symmetry breaking sets an upper bound on tan /3. 
This is mainly a consequence of two difficulties: Firstly, at large tan (3, both m 2 Hu and vn 2 Hd 
become negative around the electroweak scale and the pseudoscalar Higgs becomes tachyonic, 
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FIG. 3: Values of Ma as a function of tan/3 in the CMSSM. The SUSY parameters are fixed as 
m i/2 = 300 GeV, mo = 500 GeV and the blue lines correspond to Aq = 500 GeV blue, green lines 
to Aq = GeV and red lines to Aq = —500 GeV. Full lines correspond to \i < and dashed lines 
to fi > 0. 

as can be seen in Fig. El Secondly, the bottom Yukawa becomes non-perturbative near the 
GUT scale. Furthermore, in the region with M <C M\n there is the additional problem of 
the lighter stau becoming tachyonic. Therefore, the points of parameter space that succeed 
in generating a correct electroweak symmetry breaking with large tan (3 are strongly reduced. 
These problems are clearly softened if the entries in the down-quark Yukawa matrices are 
small as it is the case with /j > 0. In fact, for positive /x, values of tan/? up to about 60 are 
still reachable in the CMSSM. 

In the same way we must analyze the possible values of Ma- The mass of the pseudoscalar 
Higgs is approximately degenerate to the mass of the heavy scalar Higgs and both of them 
are related to the value of As already pointed out, at large values of tan/3 it is possible 
to reach smaller values of Ma due to due the fact that M\ d gets negative which also is the 
reason for the smaller value of |/i|. This can be seen in Fig. |H] where we show the possible 
values of Ma in the CMSSM as a function of tan/3 for fixed values of M1/2 = 300 GeV and 
M = 500 GeV. As expected, we see the different behaviour of M4 for positive or negative 
fi at large tan (3. For positive \x the value of Ma is larger than the corresponding value for 
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FIG. 4: Correlation between Ma and BR (B s — * All points satisfy the dark matter con- 

straint. The grey points do not survive constraints from B — > X s j and a^. Orange dots correspond 
to scenarios in which Cj S (/ib) > although they do not satisfy the BR(6 — ► sl + l~~) constraint. 

negative fi and the same value of tan f3 due to the smaller effect of the down-quark Yukawas 
on the RGE running of m^. In this plot, we do not include the experimental constraints 
on the SUSY parameter space from direct and indirect processes. In fact, making a general 
scanning of all the SUSY parameters, for fi > 0, we can still find values of Ma as low as 150 
GeV corresponding to tan (3 > 50 consistently with all constraints, see for instance Fig. 0] 

Inserting Ma = 150 GeV, tan/3 = 60 and fi > in the approximate formula given in 
Eq. (fTT]) we obtain 

A 2 (H 2 (ml /u 2 ,m 2 In 2 )) 2 
BR (B„ - u + u ) = 1.9 x 10" 4 4 , y V tJ ^,_ ^ }) ^ . (46) 



^ ^ ) ^2 (i + eo tan/?) 2 (l + e3tan/3) 2 ' 

Taking rrig ~ fi ~ m t ~ and using i7 2 (l,l) — —1/2 we have eo — 0.012 and e 3 ~ 0.015. 
Therefore, a BR(5 S — > between 10~ 7 and 5 x 10~ 6 would be the maximal value we 

can get in this model. In any case, we expect these regions of large tan/5 to be strongly 
constrained by other FCNC processes like b — > 57 and Sa^. In Figs. |U and we present 
the results for this branching ratio obtained using the complete formulae and taking into 
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FIG. 5: BR(B S -> ,ii + ju~) as a function of tan/3 in the CMSSM (upper plots) and MFV (lower 
plots). See the caption in Fig. ^ 



account all the constraints. Although many points that give rise to a large BK(B S — > /U + yU~) 
are already excluded by the stringent bounds from FCNC processes, cancelations between 
different contributions in b — > 37 still allow a large BR(5 S — > /i + /i~). 

In Figure |Hl we compare the chargino and charged Higgs contributions to the Wilson 
coefficient CV(Mjy) for the case of a BR(i? s — ► /i + /i~) just below the current experimental 
bound. With positive fi, charged Higgs and chargino contributions tend to compensate 
and we can find many points in agreement with the B — > X s 7 constraints. Moreover 
M"jj + ~ M\ + Myy, thus a small pseudoscalar mass implies also a small charged Higgs mass. 
In fact, the experimentally allowed values for C\ ot range from ~ —0.23 to ~ —0.08 while the 
contribution from the W-boson is about —0.17. This agrees with Fig. El where we can see 
that the chargino and charged Higgs contributions must be strongly correlated: in fact, their 
absolute value can be as large as 0.26 while their sum must range between -0.05 and +0.10. 
Therefore we see that, with positive fi, it is possible to have a light pseudoscalar and large 
tan/5 generating a large BK(B S — > while at the same time keep the sum of chargino 
and charged Higgs contributions to b —>■ 57 under control. Fig. El also shows a couple of 
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FIG. 6: Charged Higgs and chargino contributions to the Cj Wilson coefficient in CMSSM as 
function of tan (3 corresponding to large BR(i? s — * fj,fj,) and [i > 0. Red points refer to chargino 
and black points to charged Higgs contributions. Orange points are points where C\ ot changes sign 
although they are not allowed by the BR(6 — > sl + l~) constraint. Light grey points do not survive 
the FCNC constraints. 

points for which the chargino contribution is sufficiently large to change the whole sign of 
the total Wilson coefficient at the low scale and still remain in agreement with experiments. 
From the discussion at the end of Sec. IIII Bl it is clear that such points are excluded by the 
present measurement of BR(B — > XJ + £~) because their contributions to C 9 and C w (that 
lie within the black ellipsis drawn in Fig. EJ) are too small. 

A third process, closely related to the B s — > fifi decay, is the B s — B s mass difference. 
As discussed in section |Hj both processes involve the same FCNC neutral Higgs coupling, 
although they have a different Ma and tan/3 dependence. Note that in this case, the FCNC 
Higgs coupling enters twice as a double-penguin contribution and this implies that the change 
of sign(/i) does not affect the relative sign of the SUSY contribution to AMb s with respect 
to the SM one, that turns out to be always negative. However, the size of the FCNC Higgs 
couplings is smaller for /i > than for fx < 0. Recently the D0 collaboration was able to 
set for the first time an upper limit on AM Bs : AM Bg = (19 ± 2) ps" 1 at 90 % C.L. (31 1 . 
It is necessary to stress that this error is absolutely non- Gaussian and that there is a 5% 
probability that AMb b lies anywhere between 21 ps" 1 and infinity. Shortly afterwards the 



Unfortunately this very precise experimental measure is not fully effective in constraining 
the SUSY parameter space due to the large uncertainty in the theoretical input parameters, 



CDF collaboration improved this bound 




17.331^ (stat) ± 0.07(syst)) ps -1 . 
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FIG. 7: Correlation between ReMf 2 and BR(1? S — > See the caption in Fig. 0J 

mainly f Bs B Bs 90]. In Fig. [7|we compare Re M{ 2 with BR(5 S — > /i/z). In fact, AM Ba = 
\M* 2 \, but in our model, as in the SM, the phase of M( 2 is at the per cent level and thus it is a 
very good approximation to the total AM Bs . Moreover, plotting Re M{ 2 allows us to control 
the change of sign of this amplitude. In this plot we fix the hadronic parameters to the values 
shown in Table ITT1 As we can see in the second plot of Fig. [7| corresponding to fi > 0, in 
the CMSSM and in the parameter space allowed by FCNC constrains it is not possible 
to change significantly the value of AMg s in the region with BR(i? s — > fi + < 10~ 6 . 
Consequently no change of sign of Re M( 2 is possible in the CMSSM with \x > 0. Given the 
present experimental bound on BR(B S — > < 1CT 7 , we conclude that double penguin 

contributions to AMb s can not lower this mass difference below 17 ps _1 . To compare this 
value with the experimental measure by the CDF collaboration we must take into account 
the uncertainty in the hadronic parameters. In figure |H] we present the range of values we 
can reach varying the value of f Bs B B ^ in the interval [0.259,0.331] GeV [9J|. As we can 
see, after taking into account the theoretical uncertainty, the AM Bs constraint can never 
compete with the bound on BR(B S — > 

In the case of negative /i, we have that both eo and e,- are negative. Eqs. (JIHS1) imply 
that the diagonal Higgs couplings to the down quarks are larger than the naive expectation, 
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FIG. 8: Correlation between ReMf 2 and BR(Z? S — ► fi + fA~) including theory errors. The points 
satisfy all the constraints. The gray shaded area corresponds to theory uncertainties. The line is 
the experimental value of AMB S 



rrid/vi, and of the corresponding couplings for /i > 0. Hence, FCNC couplings are larger and 
the effects of down-quark Yukawa matrices in the RGE's is enhanced. This implies that the 
previously encountered problems with electroweak symmetry breaking occur now already at 
tan/5 ~ 45. Regarding the possible values of Ma, we see in Fig. El that, for negative /i, Ma 
tends to be much smaller than in the // > case due to the effects of the large down-quark 
Yukawas on mn d - Taking into account the constraints from BR(5 — > Xg'y) and we can 
only find allowed points with Ma > 750 GeV and tan (3 ~ 40. 

Taking // < 0, x ~ y ~ 1 and m~ g ~ ji ~ mi in Eq. we get e — —0.012 and 

63 ~ —0.015; thus, in principle, we obtain an extra enhancement from the denominator 
of Eq. (jlHl) instead of the suppression we get in the case of /i > 0. However using the 
allowed masses (as obtained above) we find that the maximum branching ratio we can get is 
about ~ 10~ 7 . In the numerical analysis, the maximum branching ratio we obtain is always 
below 10 -8 as we can see in Fig. because we have never obtained simultaneously these 
extreme values for the masses and the quantities e , Sj while being consistent with all the 
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different constraints. In particular the b — > S7 constraint requires heavy squarks and Higgs- 
bosons and Sa^ requires in additions heavy sleptons. Therefore, the SUSY contributions to 
BK(B S — > /i/i) remain small for negative /1 in the CMSSM. 

Regarding the B s — B s mass difference, the double penguin contribution has opposite 
sign to the SM contribution; hence, reducing the predicted value of AMg s . As we explained 
above, Higgs FCNC couplings are now larger for /x < and we expect larger effects than 
for /i > 0. In fact we see in Fig. [7|that with negative fi it is still possible to change the 
sign of Re M( 2 , although only very rarely. However, once we impose the indirect constraints 
we see that no change in AMb s from the SM expectations is possible. This result seems to 
differ from the analogous plots (Fig. 23) of Ref. j^]. In fact, both figures would perfectly 
agree if the b — > 57 constraints are not imposed. This is only due to the difference in the 
MSSM models considered in both works. Buras et al. analyze a general MSSM defined at the 
electroweak scale with the different parameters unrelated and only constrained by low energy 
experimental observables. In this framework they have the freedom to assume a relative sign 
between A t and /1 independently of the sign of \x. In our model, all the parameters are defined 
at Mqut and the sign of A\ at the electroweak scale is always negative due to the dominant 
(negative) gaugino contribution in the RGE evolution. Therefore our results would agree 
without the b — > 57 and <5a M constraints: taking into account these constraints, large effects 
are not possible in this GUT inspired model. 

Finally we analyze the effect of dark matter constraints on this process. Although the 
lightest neutralino of the CMSSM is a good dark matter candidate, it is well-known that in 
most of the parameter space of this model a too large dark matter density is generated. Only 
under certain specific conditions a correct value of Qy.oh 2 is generated after annihilation of 
the excess neutralinos among themselves or coannihilation with other light sparticles. The 
regions where these conditions are possible are: 

1. The bulk annihilation region, at low values of mu2 and m , where neutralinos annihi- 
late in pairs at a sufficient rate via t-channel slepton exchange. 

2. The stau coannihilation region, where neutralinos coannihilate with staus, given that 

3. The focus point region, at large m where the value of is small and the neutralinos 
have a significant higgsino component enhancing the annihilation into WW and ZZ 
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FIG. 9: Correlation between Ma and M^o. Blue and orange points corresponds to stau co- 
annihilation and funnel, respectively. In the MFV case, the dark matter constraint can be also 
fulfilled by a neutralino with a sizeable higgsino component (red points). 

pairs. 

4. The funnel region, which occurs at very large tan j3 ~ 45-60 and the pseudoscalar mass 
rriA ~ 2m^o so that neutralinos annihilate into a fermion pair through a pseudoscalar- 
Higgs boson in the s channel. 

Given that in this paper we are mainly interested in the large tan f3 region, it is clear that 
we can expect the funnel region to play a special role in our analysis. In Figs. 01 and El we 
plot the allowed points in the rriA-mx Q anc ^ m f~ m x° pl anes - Here we select points in the 
stau coannihilation region (black points) for which the difference between m? and m^o is 
less than 10%. Similarly we select the funnel points when {Ma — 2m^o) /V a < 6. Taking 
into account that Ta oc Ma tan 2 f3 {3ml + m r)' we have that at large tan/3 the difference 
between neutralino and pseudoscalar masses can be larger. In Fig. we see that the funnel 
region is wider at low pseudoscalar masses because low pseudoscalar masses correspond to 
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large tan (3. In the same way, Fig. ITUl shows clearly all the coannihilation points concentrated 
around the Ma — 2 m^o line. It is interesting to notice that we have points where both the 
funnel and coannihilation mechanisms are active (although we plot them as coannihilation 
points, i.e. black points). In fact, in these points the cross sections for x°^i ~ > 7 r+ and/or 
f+ff — > ff constitute more than 10% of the total annihilation of SUSY particles. This 
means that these points correspond to the coannihilation region. Moreover the ratio of 
X°X° ~ * bb and x°X° t~ + t~ channels is proportional to (mj/m,-) 2 which means they are 
mediated by a Higgs particle and hence they belong also to the funnel region. 

In the CMSSM we scan values of tuq and Mi/ 2 up to 1 TeV and thus we do not enter the 
focus point region. In any case, points belonging to the focus point region have very heavy 
sfermion and pseudoscalar Higgs masses and therefore they have no effect on the processes 
we are analysing here. The bulk region is already very constrained in the CMSSM scenario 
due to the bounds on the masses of the lighter chargino and the light scalar Higgs boson 
and pure bulk annihilation is nearly excluded. 

In Figure ITTI we present for completeness the correlation between Ma and M^o imposing 
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FIG. 11: Impact of the one sigma dark matter bound on the correlation between Ma and M^o. 



the full one-sigma bound on the dark matter abundance from Eq. ([44)1 . In this case the 
number of allowed points is reduced but the main features of our analysis remain valid. 



C. B s ^ fi + fi- in the general MFV MSSM 

In this section we analyze the effects of flavour changing Higgs couplings in a general MFV 
MSSM defined at the GUT scale. In this model we assign different soft masses and trilinear 
couplings to the different representations under the SM gauge group and thus we have 13 
independent parameters. The main difference with the CMSSM analyzed in the previous 
section concerns precisely the Higgs soft masses that are not related with the squark or 
slepton masses. In this way we can obtain lighter Higgs masses while at the same time the 
sfermion masses are heavy enough to satisfy the stringent FCNC constraints. 

Again we distinguish the cases of fi > and \i < 0. Analogously to the CMSSM case, 
jj, > implies that both eo and e,- are negative. Therefore Yukawa couplings are reduced 
and their effects in FCNC and the RGE evolution is smaller. 

The values that can be obtained obtain for rriA and tan f3 in this general MFV model are 
similar to the values obtained in the CMSSM case, as can be seen in Fig. Therefore, in 
principle, similar vales for the BR(5 S — > are possible in both models. However, as we 

will show here, after taking into account the FCNC constraints, combinations of (tua, tan/3) 
that would be forbidden in the CMSSM are now allowed in the MFV model. Thus it is easier 
to obtain larger values for BK(B S — ► consistently with the FCNC constraints. This 

can be seen in Figs.|Uand|3]where for positive \x we obtain similar values for BR(_B S — > fi + fi~) 
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FIG. 12: Values of uia as a function of tan/3 in the CMSSM. In this plot all sfermion masses are 
equal to 500 GeV and = 300 GeV. The red lines correspond to mn x = rriH 2 = 150 GeV, green 
lines mux = 150 and mjj 2 = 300 GeV, blue lines mu x = 300 GeV and m# 2 = 150 GeV and black 
lines to m,H\ = m H 2 = 300 GeV. Aq = GeV for all the lines. Full lines correspond to \i < and 
dashed lines to fi > 0. 

as in the CMSSM case although the dependence of BK(B S — > on Ma and tan /3 is less 

defined than in the CMSSM case. In these figures we find larger values for the branching 
ratios with respect to the CMSSM case with smaller values of Ma or tan/3. 

How the b — > 57 and a M constraints are satisfied in this case?. The main difference 
between the CMSSM and the MFV model is the fact that now different scalar masses are 
independent and therefore it is possible to have light Higgs masses with heavy sfermions. 
Apart from this difference, the case of positive \x is qualitatively similar to the CMSSM. In 
Fig. E3we present the values of the Wilson coefficients and for \x > and fi < 0. 
For n > the situation is very similar to the situation we found in the CMSSM. The only 
difference is that we can find points with BR(5 S — > > 5 x 10 -8 for comparatively 

smaller values of tan (3 and that we do find several points with C* > 1 that change the total 
sign of C 7 , although they are always forbidden by the B — > XJ + l~ constraint as expected 

EL 

Regarding the B s — B s mass difference, we see in Fig. 0that for /x > 0, the MFV case is 
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FIG. 13: Charged Higgs and chargino contributions to the CV Wilson coefficient in MFV as function 
of tan (3. The lower plots correspond to a large B s — * \x\x scenario. Red points refer to chargino 
and black points to charged Higgs contributions. Orange points are points where C\ ot changes sign 
although they are not allowed by the BR(6 — ► sl + l~) constraint. Grey points do not survive the 
FCNC constraints. 

still very similar to the CMSSM one. However we can see that now it is possible to change 
the sign of Re(M^ 2 ) if we disregard the (B s — * /i + /x~) constraint. Once we impose the new 
CDF and D0 constraint on BR(i? s — ► /x + /x~), no sizeable change in AMg s is possible and 
taking into account the theoretical errors as done in Fig. [HI the allowed points always agree 
with the recent experimental result from the CDF collaboration. 

In all these plots, the main difference between the CMSSM and the MFV cases can be 
seen for /x < 0. Once more \x < implies that both eo and 6j are negative. From Figs.|3]and|I3 
we can see that, as opposed to the CMSSM case, we can reach BR(£> S — > /x + /x_) — 10~ 6 
with /x < consistently with all the constraints and for values of tan (3 as low as 40. This 
is due to the negative sign of the parameters eo and 63. In fact for /x > the factor 
1/ ((1 + e tan (3){1 + e 3 tan (5)) can amount to a suppression by a factor 10 for tan /3 = 50 and 
eo = 0.012 while it represents an enhancement factor of 20 for tan/3 = 40 and eo = —0.012. 
Therefore, for the same values of and tan (3 we can get an enhancement of two orders 
of magnitude when changing from /x > to [L < 0. 

The main question is how are the different indirect constraints satisfied for these points. 
From Fig. El we see that, in this case, both and C7 are negative as expected and both 
of them very small as required by the b — > 57 constraint summing up to, at most, 0.05. This 
is possible only with a rather heavy spectrum although we want to keep the Higgses as light 
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as possible in order to have a large BK(B S — > Similarly, even taking a conservative 

bound at 3a, only a rather heavy SUSY spectrum is allowed. However, in this model the 
Higgs masses are independent from the rest of the SUSY spectrum and hence we can still 
find sizable Higgs contributions to the decay (B s — > fi + fi~). In any case it is clear from the 
density of point with large branching ratio that all these points require a certain degree of 
fine-tuning to satisfy these requirements. 

Finally we discuss the dark matter constraints. From Figs. 03 and El we can see the 
regions of parameter space allowed by the dark matter constraints. The main difference 
between the CMSSM case and the general MFV scenario is the appearance of a new set of 
allowed points that correspond to neutralinos with a sizeable higgsino component. This has 
also been noted previously in models where the Higgs-mass parameter has been decoupled 
from the sfermion parameters 



23, 
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] . Now the coupling of the Z-boson to an LSP pair is 
enhanced and the neutralinos annihilate through a Z-boson in the s channel. This situation 
is similar to the focus point region in the CMSSM. However, in contrast to the focus point 
situation, the sfermions and Higgs masses can be still of the order of a few hundred GeVs. 
This leads to the observed smear-out of the different regions in the MFV scenario. 



V. LARGE BK(B S -> /xV ) COLLIDER PHENOMENOLOGY 

In the first part of this section we assume that a branching ratio for the decay B s — > 
[i + fi~~ has been measured between 8 x 10~ 8 and 10 -7 and explore the consequences for the 
phenomenology of high energy colliders | . In the second part we will briefly comment on 
the occurrence of flavour violating SUSY decays. 

We have seen in the previous section, that such a large branching ratio for B s — > /i + yU~ 
requires a light pseudoscalar Higgs boson and large tan/3 as it scales as tan 6 /3/m^ - In the 
case of the CMSSM this indeed implies an upper bound of about 570 GeV and a low bound 
on tan (3 of about 50 as can be seen in the left plot of Fig. where the black points are 
those satisfying all low energy bounds while yielding the large branching ratio. In the right 
plot of Fig. El we see that these bounds get considerable weaker in the MFV case. The 
reason is that the tight correlations between the various SUSY masses is broken up and it 
is easier to satisfy the low energy constraints, in particular b — > sj. Nevertheless we find 
a lower bound of about 20 for tan/3 and an upper one of about 950 GeV for m^o. In the 
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FIG. 14: Correlation between tan/3 and Ma in a large B s — > scenario. Orange points have 
BR(B S -> n/j) < 8 x 1(T 8 ; blue points have 8 x 1(T 8 < BR(B S -> W ) < 1(T 7 . 
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FIG. 15: Correlation between and Mj± in a large -B s — ► m 1 scenario. Orange points have 
BR(B S -»/i/j)<8x 1(T 8 ; blue points have 8 x 1(T 8 < BR(B S W ) < 1(T 7 . 

CMSSM upper bounds on the masses of other particles are obtained, e.g. the chargino has 
to be lighter than 600 GeV as can be seen from Fig. where m -+ is shown versus . In 
the case of the MFV scenario the bound on m~+ is about 850 GeV. Note that the upper 
bound on m~+ implies in CMSSM (MFV) also an upper bound on m- g of about 1.7 TeV 
(2.4 TeV). As a consequence the discovery of SUSY 92j is guaranteed in both scenario due 
to gluino production. An important question is if the heavier Higgs boson can be detected 
at LHC. In the CMSSM it turns out that in the region with large BK(B S — ► /i + /i~) the 
heavier neutral Higgs bosons can be detected via their the decays into the t + t~ final state 
with 30 fb _1 , see |23j and references therein. In the region with Mao ^ 350 GeV also the 
yU + yU~ final state can be used in case of an integrated luminosity of 60 fb _1 . Also the charged 
Higgs bosons should be visible via its decay into rv T (^j]. In the case of the general MFV 
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scenarios the regions with tan/3 < 25 and/or ttlao ^ 800 GeV will most likely require larger 
statistics for detecting the heavy Higgs bosons. In the MFV scenario and to a much lesser 
in extend in the CMSSM scenario there are some points with m#+ < m t — implying 
that the branching ratio t — > br + u T will be modified. Numerically we find that there are 
3 (10) data points (out of 40000 points) in the region with BR(I? S — > > 8 • 10~ 8 ) 

(BR(B S -► < 8 • 10~ 8 )) with a maximal value of BR(t -> br + u T ) of 12 (21) %. 

We have also checked if the running from Mqut down to Mewsb can induce any sizable 
flavour violating branching ratio in the decays of SUSY particles. As expected these decay 
branching ratios are in general very small, of the order 10 -5 and below. There are a few 
points in the MFV scenario where BRs of the order 10 -3 can be obtained e.g. for b% — > Xi c 
due to a kinematical suppression of b\ — ► \[t and b\ — > xfy*- I n these points b\ is mainly 
a bi and thus the coupling b\-y§-b is reduced due to the small hypercharge of Ql implying 
a reduction in the decay mode b\ — > Xi^- However, it has to be expected that the SUSY 
background of Sl — * Xi c w ih be too large to observe this mode at LHC . 
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an approximate solution of the h-c-Xi coupling is used whi ch results 
from a one-loop integration of the flavour violating RGEs for Mq 23 and A u ^3 jlQ8| resulting 
in a mixing between cl and t\. The projection e of cl onto t\ can be written as: 

cos 9 t + sin 9 t 



2 2 
m~ — 771- 

L L 1 1 



(47) 



a 2 ( M 2 

A L = In [ -f j U;^ c6 n 2 (M| i3 + M| fl 3 + + A 2 b ) (48) 



o 2 / M 2 

^"ifc'n^j^ 2 '"" 4 ' (49) 

resulting in the ti-c-Xi interaction 

£ = -V2 + ^) ecP^it! + /i.c. (50) 

Here 9t is the mixing angle between left- and right stops. All the quantities should be taken 
at the Mqut as this is an one-step integration of the RGES. Due to the largeness of m t it 
can happen that all two-body decays of t\ are kinematically suppressed or even forbidden 
at tree-level and that the main decay mode is t% —>■ Xi c - m the MFV scenario we find 
indeed a couple of points where the tree-level decays are suppressed. These are the points 
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FIG. 16: Upper plot: approximate t\-c-Xi coupling over the exact one as a function of tan/3 in 
the CMSSM. Lower plot: exact ti-c-x? couplings a function of tan f3 in the CMSSM. The SUSY 
parameters are fixed as mj/ 2 = 300 GeV, uiq = 500 GeV and the blue lines correspond to Aq = 500 
GeV blue, green lines to Aq = GeV and red lines to Aq = —500 GeV. Full lines correspond to 
[A < and dashed lines to fj, > 0. 

shown in the lower left corner of Fig. Depending on the coupling strength it might be 
possible that t\ hadronizes before decaying and therefore it is important to check the quality 
of this approximation. In Fig. EH we show the ratio of the approximate coupling over the 
exact calculation using numerical solutions of the RGEs. One sees that for non-zero Aq 
the approximation gives roughly the order of magnitude but can be of by a factor 100 in 
certain cases. Note, that this corresponds to a factor 10 4 in the partial width. In the case 
Aq = the difference can even be larger which is due to the smallness of the exact coupling. 
The minimum of the exact coupling around tan /3 = 10 is due to a cancelation between the 
ti-CL mixing and the £r-cl which is not obtained within the approximation because the 
tR-CL is zero in this case. We have used the approximation formula for the calculation of 
the decay t\ — > c\\ getting a maximal branching ratio of 7 • 10 -3 which is about a factor 
of 200 larger than the exact calculation. From this we conclude that results obtained using 
this approximate formula have to be taken with some care. From this we conclude that an 
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observation of this decay mode t\ — > cxi would be a clear indication of the existence of 
non-minimal flavour violation structures. 

VI. CONCLUSIONS 

We have analysed flavour changing processes in a minimal flavour violating supersym- 
metric model at the GUT scale. Two examples of this kind have been studied and compared 
in more detail, the usual Constrained MSSM and a completely generic MFV model defined 
at the GUT scale. In the analysis of these models we have used two loop renormalization 
group equations taking into account for the first time the complete flavour structure in the 
running of the parameters. We have applied all the updated constraints available at the 
moment including the last measurement of the B s mass difference by CDF and the dark 
matter relic abundance from the 3rd year results of the WMAP experiment. 

The main objective of our analysis has been the study of the B s — > decay. This 

decay rate is enhanced by additional powers of tan (3 when mediated by neutral Higgses. In 
both the generic MFV and even the CMSSM scenario we indeed find that enhancements of 
the BR(5 S — > of several orders of magnitude are still possible consistently with all 

other constraints. This implies that the present bound from the D0 and CDF experiments 
can already rule out a large part of large tan (3 and small M4 region in the parameter space 
of these models. 

In the CMSSM, large values (close to the current bound) for the B s — > decay would 

necessarily imply a value of tan (3 > 50, a light pseudoscalar, M A < 570 GeV and positive 
sign of the /j, parameter. These points of the parameter space can still satisfy the stringent 
b — > «7 constraints due to the opposite sign of the chargino and charged-Higgs contributions. 
We have also checked explicitly that the new constraint from BK(B — > X s l + l~) eliminate 
all the points of the parameter space where a large chargino contribution is able to reverse 
the sign of the C-j Wilson coefficient. Unfortunately the new experimental results on the 
B s mass difference are not effective to constrain the parameter space of the CMSSM once 
the B s — > bounds are applied. This is due to the large theoretical uncertainty in the 

calculation of the B s mass difference. With respect to the dark matter constraints, we have 
seen that the allowed points with sizeable flavour changing effects correspond to the funnel 
and coannihilation regions. 
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In the generic MFV model many of the well-defined regions observed in the CMSSM 
scenario are smeared out due to the fact that the Higgs masses are now independent from 
the other soft-breaking parameters. Now, having 8 x 1CT 8 < BR(B S — > fi + '//") < 1 x 1CT 7 
would only require tan/5 > 20 and M\ < 950 GeV. Moreover, it is still possible to obtain 
these large branching ratios with negative values of the \x parameter at the expense of some 
fine-tuning to fulfil the b — * 57 constraint and satisfying the Sa^ bound at the 3 a level. As 
in the CMSSM case, the lack of more precise hadronic parameters in B s decays prevents 
us to use the precise experimental results to further constrain the allowed parameter space. 
Furthermore we have seen that in a generic MFV scenario at the GUT scale a new set of 
points with a correct dark matter abundance appears, corresponding to neutralinos with a 
large higgsino component. These points can still have relatively light Higgses and sfermions 
compared to the focus point region of the CMSSM, and therefore their contribution to FCNC 
processes can be sizeable. 

Finally we have also explored the phenomenology to be expected at collider experiments 
as LHC if a 8 x 10~ 8 < BR(_B S — > /x + < 1 x 10~ 7 is measured. The main feature of this 
scenario is that the SUSY spectrum is relatively light and SUSY should be discovered at 
LHC. The neutral and charged Higgses in this scenario should also be visible at the LHC. In 
some case with light m H + the branching ratio of the top quark into br + u T will be modified 
from the SM expectations due to charged Higgs mediation. We have studied the flavour 
violating decay t — > Xi c an d we have found that large branching ratios are not found in 
contrast with the results in the literature obtained using approximate formulas. 
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APPENDIX A: SFERMION MASS MATRIX 



The 6x6 sfermion mass matrices are given by 
/ M 



Mf 



where 



V 



I + (Tf - Q f sin 2 W ) cos 2(3 m| + M) Y£ f -j= Q((3) - M fl iQ{(3) 

Mj R + Q f sin 2 9 W cos 2(3 m| + Mj 



Y Af ^={l((3)-M f ^Q((3) 



\ 

(Al) 



(A2) 



e([3) = cot 0, tt([3) = sin (3 for Tf = \ 
Q(f3) = tan (3, Q((3) = cos (3 for Tf = -\ ' 
and Y Af are the trilinear matrices equal at M GUT to Y Af = YfAf. These matrices are 
diagonalized by the 6x6 unitary matrices Tj: 

diag(M /i ,...,M /6 ) = r / -M)-r} . 

The 6x3 left and right block components of the mixing matrices are defined as: 



(A3) 



■p6x6 _ Zp6x3 p6x3\ 
/ V fL fR J ■ 



(A4) 



In the flavour blind scenario, the most important off-diagonal entry in the above squared 
mass matrices is the third generation LR mixing. Below we present the analytic expressions 
of the 2x2 stop system: 



Ml 



Mf 

ILL 

e^iMl 

tLR 



e -mM? 

tLR 

Mf 

l RR 



(A5) 



where 



Mf 

III 



1 2 

m Q 3 + ( 2 ~ 3 sin2 cos 2 @ m z + m t > 



M Lr = m U3 + 3 sin2 6w cos 2 I 3 m z + m t , 
m t \A t -fi*Q((3)\ , 

ip- t = arg[A t -»*e(P)] , 



tLR 



The eigenvalues are given by 



Ha = ( M L + M U =f m LL - my + 



(A6) 

(A7) 
(A8) 
(A9) 

(A10) 
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with m| < m? 2 . We parametrize the mixing matrix TZ f so that 



TV 



t, 



V 



hi V 

where <^ is given in Eq. (|A9|) and 



e 2 Vi cos 0£ e 2 'Pi sin 0£ 
-e^tsin^- e^a^cos^ 



cos 61 



A 2 = (M 



-Ml 



sin ^ 



M 



III 



> 



(All) 



(A12) 



APPENDIX B: CHARGINO MASS MATRIX 



MS 



The chargino mass matrix 

M 2 mwV^ sin P 
^m w \/2cos (3 fi 
can be diagonalized by the biunitary transformation 

where U and V are unitary matrices such that are positive and < m 



(Bl) 



(B2) 



*2 



APPENDIX C: NEUTRALINO MASS MATRIX 



We define N a j as the unitary matrix which makes the complex symmetric neutralino mass 
matrix diagonal with positive diagonal elements: 



where m^o < m^o for j < k. In the basis: 

ifj a = {B,W°,H a ,H b } , 

the complex symmetric neutralino mass matrix has the form 

Mi —m z cos (5 sin Q w m z sin (3 sin 6 W 

M 2 cos /3 cos 6 W —m z sin /3 cos 6w 

—mz cos /3 sin 6 W m z cos /3 cos 9w — // 

i m z sin /3 sin 6^ —m z sin /3 cos — ji 



(CI) 



(C2) 



(C3) 
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APPENDIX D: LOOP FUNCTIONS 



In this appendix, we collect the different loop functions in the text. 



H 2 {x,y) 
S (x) 

f2(x) 

fs(x) 
Y(x) 
W(x) 

h(x) 
fe(x) 
fr(x) 

I 2 2 2\ 

c (mi,m 2 ,m 3 ) 
c 2 (ml,ml,ml) 



d 2 (ml,ml,rnl,ml) 



x logx 



+ 



2/ logy 



(l-x)(x-y) (l-y)(y-x) 
Ax — llx 2 + x 3 3a; 3 

4(1 - x) 2 2(1 -x) 3 ° g% 

-7 + 5x + 8x 2 2x - 3x 2 



6(1 -xf (1-x) 4 
3x — 5x 2 2x — 3x 2 
2(1 -xf + (l-xf ° gX 



logx 



X 



2(x-lf 
-3x 3 + 2x 2 



Ax + 3 + 2\nx) 



2{x - l) 4 
-32x A + 38x 3 



log a: + 



ix + 5x — 7x 



12(x- l) 3 
15x 2 — 18x 



18(x - l) 4 
-18x 4 + 163x 3 - 259x 2 + 108x 
36(x- l) 3 

X X 

i + 71 ^ lo g x 

1 — x (1 — xy 

38x - 79x 2 + 47x 3 

6(1 -x) 3 
52 - lOlx + 43x 2 

6(1 -x) 3 

mf log -4 



logx + 



4x — 6x 2 + 3x 4 
+ (1-x) 4 lQgX 
6 - 9x + 2x 3 



+ 



(m 2 — m 2 .)^ 2 — m|) 

A 

,2 



(I-*) 4 l0gX 
+ (mi <-> 7712) + (mi <-> 7723) 



m| log ^ 



(ml — m^mf — m§) 
4 1 m f 

™1 Wg -TjT 



(m 2 — m^mj — mj)(m 2 — m 2 ) 
+ (mi <-> 7B4) 



+ (mi <-> 7712) + (mi «-> 7713) 
+ (mi <->• m 2 ) + (mi 



(Dl) 
(D2) 
(D3) 
(D4) 
(D5) 
(D6) 



(D7) 
(D8) 

(D9) 

(D10) 

(Dll) 

(D12) 
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(D13) 



[1] For a recent review, see 

A. Masiero and O. Vives, Ann. Rev. Nucl. Part. Sci. 51 (2001) 161 |arXiv:hep-ph/0104027| ; 

42 



A. Masiero and O. Vives, New Jour. Phys. 4 (2002) 4. 
[2] G. D'Ambrosio, G. F. Giudice, G. Isidori and A. Strumia, Nucl. Phys. B 645, 155 (2002) 

|arXiv:hep-ph/0207036| . 
[3] M. Dine, R. G. Leigh and A. Kagan, Phys. Rev. D 48 (1993) 4269 |arXiv:hep-ph7 9304299 ; 

G. G. Ross, L. Velasco-Sevilla and O. Vives, Nucl. Phys. B 692 (2004) 50 

|arXiv:hep-ph/0401064| . 

For a review and further references see: 

G. G. Ross, "Models of Fermion masses", Published in TASI 2000 ed. J. L. Rosner (World 
Scientific, New Jersey, 2001) 



[4] 
[5] 
[6] 
[7] 

[8] 

[9] 

[io; 
[ii 

[12 

[13 
[14 

[15 

[16 
[17 



K. S. Babu and C. F. Kolda, Phys. Rev. Lett. 84 (2000) 228 |arXiv:hep-p h/9909476 . 

S. R. Choudhury and N. Gaur, Phys. Lett. B 451 (1999) 86 [arXiv:hep-ph/9810307| . 

G. Isidori and A. Retico, JHEP 0111 (2001) 001 |arXiv:hep-ph/0110121| . 

C. Bobeth, T. Ewerth, F. Kruger and J. Urban, Phys. Rev. D 66 (2002) 074021 

|arXiv:hep-ph/0204225| . 

P. H. Chankowski and L. Slawianowska, Phys. Rev. D 63 (2001) 054012 
|arXiv:hep-ph/0008046| . 

A. J. Buras, P. H. Chankowski, J. Rosiek and L. Slawianowska, Phys. Lett. B 546 (2002) 96 
|arXiv:hep-ph/020724l'| . 

A. J. Buras, P. H. Chankowski, J. Rosiek and L. Slawianowska, Nucl. Phys. B 659 (2003) 3 
|arXiv:hep-ph/0210145| . 



G. Isidori and P. Paradisi, arXiv:hep-ph/0605012 



M. Carena, A. Menon, R. Noriega-Papaqui, A. Szynkman and C. E. M. Wagner, 
|arXiv:hep-ph/0603106[ 

S. Baek, P. Ko and W. Y. Song, JHEP 0303 (2003) 054 |arXiv:hep-ph/0208112| . 

A. Dedes, H. K. Dreiner and U. Nierste, Phys. Rev. Lett. 87 (2001) 251804 

|arXiv:hep-ph/0108037| . 

R. Arnowitt, B. Dutta, T. Kamon and M. Tanaka, Phys. Lett. B 538 (2002) 121 
|arXiv:hep-ph/0203069| . 

T. Ibrahim and P. Nath, Phys. Rev. D 67 (2003) 016005 |arXiv:hep-p h/0208142 . 

J. K. Mizukoshi, X. Tata and Y. Wang, Phys. Rev. D 66 (2002) 115003 



arXiv:hep-ph /0208078 . 

43 



A. Dedes, H. K. Dreiner, U. Nierste and P. Richardson, arXiv:hep-ph/0207026 

H. Baer, C. Balazs, A. Belyaev, J. K. Mizukoshi, X. Tata and Y. Wang, JHEP 0207 (2002) 

050 | arXiv:hep-ph/0205325| . 

T. Blazek, S. F. King and J. K. Parry, Phys. Lett. B 589 (2004) 39 |arXiv:hep-ph7 0308068 . 
G. L. Kane, C. Kolda and J. E. Lennon, |arXiv:hep-ph/03 10042 



J. R. Ellis, K. A. Olive and V. C. Spanos, Phys. Lett. B 624 (2005) 47 
|arXiv:hep-ph/0504196| . 

S. Baek, Y. G. Kim and P. Ko, JHEP 0502 (2005) 067 |arXiv:hep-ph/0406033| . 

R. Dermisek, S. Raby, L. Roszkowski and R. Ruiz De Austri, JHEP 0304 (2003) 037 

|arXiv:hep-ph/0304101| . 

J. Ellis, K. A. Olive, Y. Santoso and V. C. Spanos, |arXiv:hep-ph/0603136 



A. Bartl, T. Gajdosik, E. Lunghi, A. Masiero, W. Porod, H. Stremnitzer and O. Vives, Phys. 
Rev. D 64, 076009 (2001) |arXiv:hep-ph/0103324| . 

N. Polonsky and A. Pomarol, Phys. Rev. D 51, 6532 (1995) |arXiv:hep-ph/9410231| . 

L. E. Ibahez, D. Lust and G. G. Ross, Phys. Lett. B 272 (1991) 251 |arXiv:hep-th/9109053| ; 

G. G. Ross, |arXiv:hep-ph/0411057l 



S. Martin and M. Vaughn, Phys. Rev. D50, 2282 (1994); Y. Yamada, Phys. Rev. D 50, 3537 

(1994); I. Jack, D.R.T. Jones, Phys. Lett. B333 (1994) 372. 

D. N. Spergel et al, arXiv:astro- ph/0603449 . 

V. Abazov [DO Collaboration] , |arXiv:hep-ex/0603029l 

" Measurement of the Bs - Bs Oscillation Frequency" , available at CDF's B-physics webpage: 
http:/ /www-cdf. fnal.gov/physics/new/bottom/bottom.html 



G. W. Bennett et al. [Muon g-2 Collaboration], Phys. Rev. Lett. 92, 161802 (2004) 
|arXiv:hep-ex/0401008| . 

M. Iwasaki et al. [Belle Collaboration], arXiv:hep-ex/0503044 



B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 93 (2004) 081802 



arXiv:hep-ex/0404006 . 

S. Eidelman et al. [Particle Data Group], Phys. Lett. B 592, 1 (2004). 
[CDF Collaboration], |hep-ex/0507091| 
CDF Collaboration, Public note 8176. 
DO Collaboration, note 5009-CONF. 



44 



[39] A. Dedes and A. Pilaftsis, Phys. Rev. D 67 (2003) 015012 |arXiv:hep -ph/0209306 . 

[40] A. Brignole and A. Rossi, Phys. Lett. B 566 (2003) 217 |arXiv:hep-ph/0304081| . 

[41] P. Paradisi, JHEP 0602 (2006) 050 |arXiv:hep-ph/0508054| . 

[42] A. Masiero, P. Paradisi and R. Petronzio, arXiv:hep-ph/0511289 

[43] A. M. Curiel, M. J. Herrero, W. Hollik, F. Merz and S. Penaranda, Phys. Rev. D 69 (2004) 



075009 |arXiv:hep-ph/ 0312135 . 
[44] S. Kanemura, K. Matsuda, T. Ota, T. Shindou, E. Takasugi and K. Tsumura, Phys. Lett. B 

599 (2004) 83 |arXiv:hep-ph/0406316| . 
[45] T. Ibrahim, P. Nath and A. Psinas, Phys. Rev. D 70 (2004) 035006 |arXiv:hep-ph/0404275| . 
[46] T. Ibrahim and P. Nath, Phys. Rev. D 71 (2005) 055007 |arXiv:hep-ph/0411272| . 



[47] W. Hollik, S. Penaranda and M. Vogt, arXiv:hep-ph/0511021 



[48] A. Abulencia et al. [CDF Collaboration], Phys. Rev. Lett. 95 (2005) 221805 [Erratum-ibid. 



95 (2005) 249905] arXiv:hep-ex/0508036 ; 

V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 94 (2005) 071802 
|arXiv:hep-e^7041 0039 . 



[49] G. Isidori and P. Paradisi, arXiv:hep-ph/0601094 



[50] A. J. Buras, P. H. Chankowski, J. Rosiek and L. Slawianowska, Nucl. Phys. B 619 (2001) 

434 |arXiv:hep-ph/0107048| . 
[51] P. Gambino, U. Haisch and M. Misiak, Phys. Rev. Lett. 94, 061803 (2005) 

|arXiv:hep-ph/0410155| . 
[52] T. Huber, E. Lunghi, M. Misiak and D. Wyler, arXiv:hep-ph/0512066 



[53] M. Ciuchini, G. Degrassi, P. Gambino and G. F. Giudice, Nucl. Phys. B 527, 21 (1998) 
|arXiv:hep-ph/9710335| . 

[54] A. J. Buras, A. Czarnecki, M. Misiak and J. Urban, Nucl. Phys. B 631, 219 (2002) 



|arXiv:hep-ph /0203135 ; A. J. Buras and M. Misiak, Acta Phys. Polon. B 33, 2597 (2002) 

|arXiv:hep-ph/020713T| . 
[55] F. M. Borzumati and C. Greub, Phys. Rev. D 58 (1998) 074004 |arXiv:hep-ph7 9802391 . 
[56] F. M. Borzumati and C. Greub, Phys. Rev. D 59 (1999) 057501 |arXiv:hep-ph/9809438| . 
[57] C. Bobeth, M. Misiak and J. Urban, Nucl. Phys. B 567 (2000) 153 |arXiv:hep-ph/9904413| . 
[58] G. Degrassi, P. Gambino and G. F. Giudice, JHEP 0012, 009 (2000) |arXiv:hep-ph/0009337| . 
[59] M. Carena, D. Garcia, U. Nierste and C. E. M. Wagner, Phys. Lett. B 499, 141 (2001) 



45 



|arXiv:hep-ph /0010003 . 
[60] G. Degrassi, P. Gambino and P. Slavich, arXiv:hep-ph/0601135 
[61] CLEO Collaboration (S. Chen et al), Phys. Rev. Lett. 87, 251807 (2001). 
[62] Belle Collaboration (K. Abe et al), Phys. Lett/ B 511, 151 (2001). 
[63] Belle Collaboration (P. Koppunberg et al), Phys. Rev. Lett. 93, 061803 (2004). 
[64] BaBar Collaboration (B. Aubert et al), Phys. Rev. D 72, 052004 (2005). 



[65] BaBar Collaboration (B. Aubert et al), hep-ex/0507001 



[66] Heavy Flavor Averaging Group, arXiv:hep-ex/0603003 



[67] T. Hurth, E. Lunghi and W. Porod, Nucl. Phys. B 704, 56 (2005) |arXiv:hep-ph7031 2260 . 
[68] P. L. Cho, M. Misiak and D. Wyler, Phys. Rev. D 54, 3329 (1996) |arXiv:hep-ph/9601360| . 
[69] S. Bertolini, F. Borzumati, A. Masiero and G. Ridolfi, Nucl. Phys. B 353, 591 (1991). 
[70] H. E. Haber and G. L. Kane, Phys. Rept. 117, 75 (1985). 

[71] T. Moroi, Phys. Rev. D 53, 6565 (1996) [Erratum-ibid. D 56, 4424 (1997)] 

|arXiv:hep-ph/9512396l . 
[72] G. W. Bennett [Muon Collaboration], |arXiv:hep^ex7 0602035 
[73] F. Jegerlehner, Nucl. Phys. Proc. Suppl. 126 (2004) 325 |arXiv:hep-ph/0310234| . 
[74] [ALEPH Collaboration], |arXiv:hep-ex/0602042| 

[75] J. R. Ellis, K. A. Olive, Y. Santoso and V. C. Spanos, Phys. Lett. B 565, 176 (2003) 
|arXiv:hep-ph/0303043l . 

[76] G. Belanger, S. Kraml and A. Pukhov, Phys. Rev. D 72, 015003 (2005) 



|arXiv:hep-ph /0502079 . 
[77] P. B. Mackenzie, |arXiv:hep^ ph/0606034 
[78] D. M. Pierce et al., Nucl. Phys. B 491 (1997) 3. 
[79] W. Porod, in preparation 

[80] G. Degrassi, P. Slavich and F. Zwirner, Nucl. Phys. B 611 (2001) 403. 

[81] A. Brignole, G. Degrassi, P. Slavich and F. Zwirner, Nucl. Phys. B 631 (2002) 195. 

[82] A. Brignole, G. Degrassi, P. Slavich and F. Zwirner, Nucl. Phys. B 643 (2002) 79. 

[83] A. Dedes, G. Degrassi and P. Slavich, Nucl. Phys. B 672 (2003) 144. 

[84] A. Dedes and P. Slavich, Nucl. Phys. B 657 (2003) 333. 

[85] B. C. Allanach, A. Djouadi, J. L. Kneur, W. Porod and P. Slavich, JHEP 0409 (2004) 044. 
[86] C. Bobeth, A. J. Buras, F. Kruger and J. Urban, Nucl. Phys. B 630 (2002) 87 



46 



|arXiv:hep-ph/0112305| . 

T. Goto, Y. Okada, Y. Shimizu and M. Tanaka, Phys. Rev. D 55 (1997) 4273 [Erratum-ibid. 
D 66 (2002) 019901] |arXiv:hep-ph/9609512| . 

S. Baek, T. Goto, Y. Okada and K. i. Okumura, Phys. Rev. D 64 (2001) 095001 
|arXiv:hep-ph/0104146l . 

T. Ibrahim and P. Nath, Phys. Rev. D 62 (2000) 015004 |arXiv:hep- ph/9908443 . 
P. Ball and R. Fleischer, |arXiv:hep-p"h"70 604249 

A. Dedes and B. T. Huffman, Phys. Lett. B 600 (2004) 261 arXiv:hep-ph/0407285| . 
S. Asai [ATLAS and CMS Collaborations], Eur. Phys. J. directC 4S1 (2002) 17. 
A. Djouadi, |arXiv:hep-ph/0503173 



H. Baer, J. Sender and X. Tata, Phys. Rev. D 50 (1994) 4517 |arXiv:hep-ph7 9404342 . 
T. Kon and T. Nonaka, Phys. Rev. D 50 (1994) 6005 |arXiv:hep-ph/9405327| . 
J. Sender, Phys. Rev. D 54 (1996) 3271. 

W. Porod and T. Wohrmann, Phys. Rev. D 55 (1997) 2907 [Erratum-ibid. D 67 (2003) 



059902] arXiv:hep-ph/9608472 . 

M. Hosch, R. J. Oakes, K. Whisnant, J. M. Yang, B. 1. Young and X. Zhang, Phys. Rev. D 
58 (1998) 034002 |arXiv:hep-ph/9711234| . 

W. Porod, Phys. Rev. D 59 (1999) 095009 |arXiv:hep-ph/9812230| . 

C. Boehm, A. Djouadi and Y. Mambrini, Phys. Rev. D 61 (2000) 095006 
|arXiv:hep-ph/9907428| . 

D. Restrepo, W. Porod and J. W. F. Valle, Phys. Rev. D 64 (2001) 055011 
|arXiv:hep-ph/0104040l . 

A. Djouadi, M. Guchait and Y. Mambrini, Phys. Rev. D 64 (2001) 095014 
|arXiv:hep-ph/0105108| . 

S. P. Das, A. Datta and M. Guchait, Phys. Rev. D 65 (2002) 095006 |arXiv:hep-ph/0112182| . 
W. Porod, Comput. Phys. Commun. 153 (2003) 275 |arXiv:hep-ph/0301101| . 



S. P. Das, A. Datta and M. Guchait, Phys. Rev. D 70 (2004) 015009 arXiv:hep-ph/0309168 . 
M. Muhlleitner, A. Djouadi and Y. Mambrini, Comput. Phys. Commun. 168 (2005) 46 
|arXiv:hep-ph/031 11671 . 

[107] T. Han, K. I. Hikasa, J. M. Yang and X. m. Zhang, Phys. Rev. D 70 (2004) 055001 
|arXiv:hep-ph/0312"T29"l . 



47 



[108] K. Hikasa and M. Kobayashi, Phys. Rev. D 36 (1987) 724. 



48 



